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The vacuum of Yang-Mills theories can be imagined as a magnetically disordered medium with 
domain structure, with color magnetic flux in each domain quantized in units corresponding to the 
gauge group center. This model leads to the prediction of Casimir scaling, i.e. the proportionality 
of string tensions of potentials (at intermediate distances) between color sources from higher- 
representations to eigenvalues of the quadratic Casimir operator. I present evidence for Casimir 
scaling in G2 lattice gauge theory. I also discuss support for some ingredients of the model from 
the recently conjectured form of the Yang-Mills ground-state wave-functional. 
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1. Introduction: the Casimir scaling hypothesis 

The explanation of how quark and gluon confinement arises in non-Abelian gauge theories 
remains a challenging problem. To make the problem simpler, one can leave out dynamical quark 
degrees of freedom and ask e.g. about the behavior of potentials between static color sources from 
different representations of the gauge group, but even this question does not have a fully satisfactory 
answer. The potentials are expected to behave differently in three regions of quark-antiquark dis- 
tances (the situation is sketched in Fig. [I] for the case of SU(2)). At short distances the interaction is 
governed by perturbative one-gluon exchange and the potentials are Coulomb-like. At asymptotic 
distances, the color charges of higher-representation static sources can be screened by an appropri- 
ate number of gluons; in SU(A0 gauge theory, potentials of all zero Af-ality representations become 
asymptotically flat, while asymptotic string tensions for representations with nonzero Af-ality will 
be the same as for the lowest representation with the same Af-ality. Most interesting is the region 
of intermediate distances, from the onset of confinement to the onset of screening. There string 
tensions depend on color-charge representations in a more intricate way: that dependence carries 
an imprint of the underlying mechanism of confinement. 

It was argued that intermediate string tensions of potentials between higher-representation 
quarks and antiquarks should be proportional to the quadratic Casimir. This Casimir scaling hy- 
pothesis ^ can be supported e.g. by large-Af factorization [§] or dimensional reduction arguments 
[@— g]. Regardless of whether one is convinced by these arguments or not, Casimir scaling is an 
experimental fact: it was observed convincingly in numerical simulations of both SU(2) and SU(3) 
lattice gauge theories (most recently in Refs. [|6]-||]). 

Any viable model of the confinement mechanism should be able to explain both Casimir scal- 
ing at intermediate and Af-ality dependence at large distances. It is not easy at all to understand 
these effects in terms of vacuum fluctuations which dominate the QCD functional integral [Qj. 

In the present contribution I will first discuss a simple model of how (approximate) Casimir 
scaling and Af-ality dependence arise within the center- vortex picture of color confinement, and then 
extend the model to non-Abelian theories with trivial center (Section |||). Its specific prediction is 
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Figure 1: Sketch of the expected behavior of static potentials in SU(2) gauge theory. 
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Casimir scaling for potentials in (center-trivial) G2 gauge theory. I will report results of numerical 
simulations that confirm the prediction (Section |||). Finally, I will outline a recent suggestion for 
an approximate wave-functional of the Yang-Mills theory that incorporates (some) elements of the 
proposed vacuum model (Section ||). Section || summarizes conclusions. 



2. A simple model of Casimir scaling and color screening 

iV-ality (or "representation class" in mathematical language) is related to the transformation 
properties of SU(/V)-group representations with respect to the group center, Zjy. The Af-ality de- 
pendence of string tensions of different-representation potentials at asymptotic distances can be 
easily understood in terms of the picture of the QCD vacuum as a center- vortex condensate, where 
the dominant field configurations are directly related to elements of the group center [jlO|]. The 
asymptotic string tension o r can be extracted from the area-law fall-off of Wilson loops in the r-th 
representation of the color group. I will illustrate the idea on the SU(2) group; its representations 
are labeled by a "spin" index j with half-integer or integer values. Let us assume the vacuum is 
filled with percolating thin center vortices and divide the loop ^ into small patches (with unit area 
for simplicity). Let p denote the probability that a patch is pierced by a vortex and assume piercings 
are random and uncorrected. It is then a simple exercise to show that: 

Wj(V)~ [p-(-\) 2j +{\-p)-{+l)Y*^\ stiff) ... minimal area of the loop <af, (2.1) 

aJ __ h[I _, + ( _ 1) «,] M /*---l-"*r'A (2.2) 
L J ... integer j. 

The asymptotic string tensions are zero for all integer-7 representations (i.e. those with iV-ality, or 
"biality", equal to 0), and nonzero and equal for all half-integer j (with iV-ality equal to 1). The 
argument can simply be extended to arbitrary N >2. 

An explanation of Casimir scaling at intermediate distances within the center-vortex model 
was suggested by Faber et al. [II]. 1 It is attributed to the fact that center vortices in the QCD 
vacuum do have finite, relatively large thickness. 2 If the cross-section of the vortex with the plane 
of the Wilson loop is fully outside the loop, the vortex does not influence its value. If it is fully 
inside the loop, the value is multiplied by (— 1) 2; in SU(2). The effect of a vortex whose core is 
only partially contained in the area of the loop is modeled by multiplication of the loop by a group 
element (in the j-th representation) that depends on a certain angle a and has random orientation 
in color space. The angle interpolates between (if the core lies fully outside the loop) and 271 (if 



it lies fully inside). This very simple model predicts, as above, ( |2.2[ ) for asymptotic distances, and 
(to a good degree of approximation) 

<3j ~ j(j +1) • • • Casimir scaling (2.3) 

for intermediate distances. The model can be generalized to N > 2, at the price of introducing more 
parameters (angles). 3 



1 A related proposal was advanced long ago by Cornwall in a proceedings contribution jl^ ] . 

independent ways of estimating the thickness of center vortices are discussed e.g. in Section 4 of Ref. piy. 

3 See Ref. p4[ ] for particular proposals for SU(3) and SU(4). 
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Recently, much attention has been paid to the Yang-Mills theory with the G2 gauge group [|T5|] . 
The exceptional Lie group G2 has a trivial center; the center of its universal covering group is also 
trivial, and there seems to be no reason to believe that the group center plays any important role 
in G2 YM theory at all. The theory is at most temporarily confining, since potentials for static 
color charges from any representation, including the fundamental one, must be asymptotically 
flat. The reason is that any color charge, even a "quark" from the (7-dimensional) fundamental 
representation of G2 can be screened by a bunch of "gluons" (from the 14-dimensional adjoint 
representation). This fact does not really contradict the center-vortex confinement scenario: G2 
does not possess nontrivial center vortices, and therefore the asymptotic string tension is zero. But 
even in the G2 gauge theory, one expects that the static potentials do grow linearly over a certain 
range of distances, from the scale where perturbation theory breaks to the onset of screening. The 
linear rise of the fundamental potential in this intermediate region was demonstrated in numerical 



simulations fll6Q. 

Dimensional reduction can be invoked to argue that even in G2 YM theory one would expect 
Casimir scaling of potentials at intermediate distances. But if, in the language of field configura- 
tions, finite thickness of center vortices was responsible for Casimir scaling in SU(A0, what should 
it be attributed to in a theory with trivial center? A common origin of Casimir scaling for "center- 



poor" and "center-rich" gauge models was suggested in Ref. [16]. We extended and improved the 



model of Ref. [11]. We assume the Yang-Mills vacuum state has a domain structure, with the color 
magnetic flux in each domain quantized in units of the gauge group center (be it trivial or not), and 
Casimir scaling results from random spatial fluctuations of the flux in each domain. 

In this model, it is assumed that if we take a 2D slice of the four-dimensional volume, we 
can split it into domains ("patches") of a typical area A^. Within each domain color magnetic 
fields fluctuate randomly and (almost) independently, with a short length of correlation /. Each 
domain contains small independently fluctuating subregions of area I 2 <C A^. The only constraint 
on fluctuating fields bounds the total magnetic flux over each domain to correspond to an element 
of the gauge-group center. In SU(2), there are domains of the center-vortex type and of the vacuum 
type. The former correspond to the nontrivial center element, —I, and represent a cross section of 
a thick center vortex; the latter carry a zero total magnetic flux. In G2 all domains will be of the 
vacuum type, since the center contains the identity element only. 

I will only summarize conclusions of the model: If the center of the gauge group contains N 
elements, there are ,/V types of domains enumerated by the value k £ {0,1,.. . ,N — 1}, each assumed 
to appear with a probability centered at any given plaquette in the plane of the loop. The effect 
of a domain (a 2D cross section of the kth vortex) on a planar Wilson loop is to multiply the loop 
by a group element 



&(a (k) ,y) =y exp 



W-jrl y\ (2.4) 



where {Mf\ are generators of the Cartan subalgebra, 5? is a random element of the group, and 
angles a'*' depend on the location of the vortex/domain with respect to the loop. If the domain is 
all contained within the loop area, then 

exp[/a w -jrl =zkl, (2.5) 
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where zk is the kth center element, I is the unit element. If the domain is outside the loop, it has no 
effect, i.e. 

I. (2.6) 



exp 



For a Wilson loop from the representation r, the averaged contribution of a domain is 

ld r , (2-7) 



where d r is the dimension of the representation r and is the d r x d r unit matrix. 

One further assumes that probabilities to find domains of any type centered at two different 
plaquettes are independent, and that for loops smaller than the typical size of the domain, the r.m.s. 
of phases a is proportional to the area of the vortex contained in the interior of the loop. Then, 
both in SU(A0 and G2, the static potential V r (R) of the representation r will be linearly rising for 
distances / <C R <C y/Ad, with a string tension approximately proportional to its quadratic Casimir: 

<3 r ^C r ... intermediate distances. (2.8) 

For very large Wilson loops most vortices will be contained within the loop, the average phases are 
proportional to the corresponding total magnetic flux through the domain, and the prediction is, for 

o r ~ F[k r , {fk}] ■■■ asymptotic distances (2.9) 
for SU(/V), where k r is the Af-ality of the representation r, and 

a r = ... asymptotic distances (2.10) 

for the G2 gauge group. 

3. Casimir scaling in G2 gauge theory 

The described domain model of the YM vacuum with magnetic disorder, containing only a 
few adjustable parameters, predicts Casimir scaling of higher representation potentials at interme- 
diate distances not only for SU(/V), but also for G2 gauge theory. The prediction can be verified 
or disproved in numerical simulations of the theory in lattice formulation. This is in principle a 
straightforward task, but not cheap for G2: simulations are quite demanding on computer resources, 
determination of string tensions requires all machinery developed in the past for calculating poten- 
tials (anisotropic lattices, ground-state overlap enhancement, smearing) plus some bits and pieces 
of information from group theory. Technical details of our calculations can be found in Ref. |/Jl, 
I will present here only a subset of representative results. 

We simulated G2 lattice gauge theory with the Wilson action on anisotropic L 3 x (2L) lattice: 

5 = -§UolReTr^oW]-^- £ Re Tr 1 (3.1) 

' [ x,i>0 W V.I /' J 

with the bare-anisotropy parameter £0 tuned so that the physical anisotropy £ = a^ hys /aj* ys equals 2. 
To increase overlap of the trial quark-antiquark state with the ground state, we used in the computa- 
tion of Wilson loops spatial links smeared by a G2 generalization of the stout smearing method [|T8|]. 
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Most of our results come from 14 3 x 28 lattice at three values of the coupling j3 = 9.5,9.6,9.7 on 
the weak-coupling side of the crossover region observed in Ref. [fill]. 

The potential between a static quark and antiquark from the representation {D} of the G2 
group can be determined from values of Wilson loops W^ D y(r,t) by a fit of the form 

- In W {D} (r,t)= C {D} + V {D} (r) • t (3.2) 

in an appropriate interval (t m i n ,t max ). The resulting potentials are then parametrized by the usual 
Coulomb plus linear form: 

V {D} (r) = c {D} -°^ + a {D} r. (3.3) 

G2 irreducible representations {D} are labeled by two Dynkin coefficients [Ai , A2] , the dimen- 
sion of the representation is given by JI^ ] : 

D = ^$-%){%-%)(%-%), (3-4) 

where l\ = 5(1 + M), ^2 = 5 (4 + X\ +3A2), and £3 = |(5 + 2Ai +3A2). The ratio of eigenvalues 
of the quadratic Casimir operator is JI^]: 

^} = ^ = J(^+4+4-y), (3-5) 

where Cp is the quadratic Casimir for the fundamental representation. 

Results for string-tension ratios for color sources from five G2 representations are summarized 
in Table [lj potentials are displayed in Fig. || for j8 =9.6 and 9.7. Fundamental and adjoint po- 
tentials for all three couplings j3 are shown together in Fig. [3| expressed in physical units defined 
by the fundamental string tension. Our results convincingly demonstrate (approximate) Casimir 
scaling for static potentials between color charges from various representations of G2. 4 The agree- 
ment between measured values of intermediate string tensions with predictions based on values of 
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Figure 2: [Left:] Potentials for different representations vs. the dimensionless r; 14 3 x 28 lattice, j3 = 9.6. 
[Right:] Ratios V^oj/Vp for different representations {D} as functions of the dimensionless r; 14 3 x 28 
lattice, j3 = 9.7. Horizontal lines show Casimir-scaling predictions. 



4 Casimir scaling of G2 potentials has recently been seen also in numerical simulations of G2 lattice gauge theory in 
(2+1) dimensions lEOO. 
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Table 1: Ratios Os d \/(Jf for the adjoint (labeled A), {27}, {64}, {77}, and {77'} representations. The last 
line is the Casimir-scaling prediction. 
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Figure 3: The fundamental (left) and adjoint (right) representation potential in physical units, all three j3 
values, from 14 3 x 28 lattice. (The constant c, cf. Eq. (3.3), is subtracted.) 



quadratic Casimirs is quite striking; they differ by at most 10-15%, and this can hardly be just a 
numerical coincidence. The results of course cannot prove that the model described in Section || is 
right, but combined with the solid evidence for Casimir scaling in SU(2) and SU(3), they provide 
support for its main ingredients (common for all groups) - a magnetically disordered vacuum with 
a domain structure. 



4. A suggestion for an approximate vacuum wave-functional of Yang-Mills theory 
in (2 + 1) dimensions 

Can one derive (at least some) elements of the described picture from first principles? Re- 
cently, we have suggested an approximate form of the vacuum wave-functional which solves the 
SU(2) Yang-Mills Schrodinger equation in the temporal gauge [21]. 



¥o[A] = exp 



^Jd 2 xd 2 y B a (x) '< 



ab 



yJ-^-Xo + m 2 



B\y) 



(4.1) 



Here B"(x) = F" 2 (x) is the color magnetic field strength, the covariant laplacian in the adjoint 
representation: 

2 



k=\ 



(4.2) 



5 Expressions below are all assumed to be properly denned on a lattice, with lattice spacing serving as regulator, but 
for simplicity I will mostly use continuum notations. 
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where ^(x) are the link fields in the adjoint representation, Xq denotes the lowest eigenvalue of 
(— 3> 2 ), and m is a constant (mass) proportional to g 2 ~ 1//3. 



The proposed wave-functional ( J4. 1[ ) is reminiscent of Samuel's Q22|], the difference is that in 
his proposal a single free parameter hiq replaces our (— Ao + m 2 ) in the denominator. The reason 
for subtracting the lowest eigenvalue from the operator (— S> 2 ) is that our numerical simulations 
indicate that its spectrum might diverge in the continuum limit. 

This wave-functional has quite a few attractive properties: 

1. In the free-field limit (g — ► 0), the covariant laplacian turns into ordinary laplacian, Xq and m 
go to 0, and *Fo[Al becomes the well-known vacuum wave-functional of electrodynamics: 

y [A] = expl-yd 2 xd 2 y[VxA^x)}^^=^ [VxA b (y)]\, (4.3) 



2. 



5. 



Eq. ( 4. 1 ) is a good approximation to the uue vacuum also in completely different corner of 
the configuration space, namely if we restrict to fields constant in space and varying only in 
time. In the large-volume limit the solution of the YM Schrodinger equation in that case is, 
up to l/V corrections: 



y Q = exp[-VR ] =exp 



1 y (A\ xA 2 )-(A l xA 2 ) 

2 y/Ai ■A l +A 2 -A 2 



(4.4) 



and exactly the same expression follows from (4T) assuming |gAi,2| S> m, y/J®. 



3. If we divide the field strength B"(x) into "fast" and "slow" components, the portion of the 
(squared) vacuum wave-functional in B s iow is approximately 



: exp 



1 



d 2 xB« low (x) B a slow (x) 



(4.5) 



This is exactly the probability measure for YM theory in two Euclidean dimensions, which 
(i) is confining for m > 0, and (ii) exhibits Casimir scaling for string tensions of all color- 
charge representations. The fundamental string tension is easily computed as Cf{P) = 
3m/ (4/3). The last expression can be used to fix the value of m in Eq. ( |4.1| ) at a given j8 
from the known value of Op. 



Confinement requires m to be larger than zero; if one takes m in the wave-functional ( |4.1[ ) as a 
variational parameter and computes (approximately) the expectation value of the Yang-Mills 
Hamiltonian, one finds that non-zero (and finite) value of m is energetically preferred. 

If we fix the mass m in the wave-functional to get the right string tension Op at a given j8, 
we can test our proposal by calculating e.g. the mass gap of the theory. We have proposed 
a recursive procedure for generating independent lattice configurations with the probability 



distribution given by the square of the wave-functional (4.1), the interested reader should 



consult Ref. [21] for its detailed exposition. We call two-dimensional lattice configurations 
obtained in this way "recursion lattices". One can compute various observables with these 
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Figure 4: Mass gaps extracted from recursion lattices at various lattice couplings, compared to + glueball 
masses in 2 dimensions obtained in Ref. [ [24] ] (denoted "expt") via standard lattice Monte Carlo methods. 



lattices, and compare results with "MC lattices", i.e. two-dimensional slices of lattices gen- 
erated in a full D = 3 lattice Monte Carlo simulation. To extract the mass gap, one computes 
the correlator 

9(x-y) = ((B"B a ) x (B h B b ) y ) - {{B a B a ) x )\ (4.6) 
and determines the mass gap M by fitting the result to the form 

7 e- MR 

(const.) (1 + iM/?) 2 -^, R=\x-y\. (4.7) 

Fig. |] shows our values for the mass gap versus /3 extracted from recursion-lattice data. It 
turns out that, given the asymptotic string tension as input, we can compute the mass gap 
fairly accurately from our wave-functional. The discrepancy from Monte Carlo results of 
Meyer and Teper [Q for the + glueball mass are at the level of at most a few (< 6) percent. 



6. The YM wave-functional in Coulomb gauge can be obtained by restricting the temporal- 
gauge wave-functional to transverse gauge fields. We have computed the color-Coulomb 
potential and the Coulomb-gauge ghost propagator using recursion and MC lattices and the 



results are in reasonable agreement, see Jeff Greensite's talk at this workshop [23]. 



The dimensional reduction form ( |4.5| ) at large distances implies an area law fall-off for large 
Wilson loops, and also Casimir scaling of higher-representation Wilson loops. The question is how 
Casimir scaling turns into AT-ality dependence, how color screening enters the game within this 
setting. It may be necessary to introduce into our wave-functional additional term(s), e.g. a gauge- 
invariant gluon-mass term advocated by Cornwall p5j]. However, there are indications that terms 



needed for color screening might be contained in our simple wave-functional and would appear as 



corrections to the dimensional -reduction form (4.5) 



Let us write the lattice vacuum state ^[i/] in the form exp (R[U]). Greensite [^] developed 
a strong-coupling technique for calculating R[U) in Hamiltonian lattice gauge theory. R[U) is ex- 
pressed as a sum over spacelike Wilson loops and products of loops on the lattice, schematically: 
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(4.8) 

The first few coefficients c,- for SU(2) lattice gauge theory in (2 + 1) dimensions were computed in 



Ref. Q26| ] and for smoothly varying fields one gets 

2 



¥<,[!/] =exp 



J d 2 x (aK B 2 -a 3 K 2 B(-£> 2 )B+...) 



(4.9) 



where a is the lattice spacing, Kq = ^cq + 2(c\ + c 2 + C3), Ki = \c\,cq'v& @{fi 2 ), and ci,c 2 ,c^ are 
<^(/3 4 ). The leading correction to dimensional-reduction form is contained in the term proportional 
to K 2 , and comes from the 1x2 loop in R[U] proportional to c\. If we evaluate the adjoint Wil- 
son loop WaC^) in the strong-coupling expansion, the c\ term provides the leading perimeter-law 
contribution, sketched in Fig. ||, proportional to (ci/2)*^*^ -4 . So the term that gives the leading 
correction to the dimensional-reduction form is also responsible for screening of adjoint loops. 



Figure 5: Tiling/screening an adjoint Wilson loop by 1 x 2 rectangles. The loop is denoted by a heavy solid 
line. 



If we now work out the leading correction to the dimensional-reduction form ( |4.5| ) by expand- 

-^ 2 -Ao. 



ing the kernel in our wave-functional ( |4.1[ ), we get: 

1 







exp 



m 



d 2 x 



^slow^slow " 



B 



slow" 



2m 2 



"^slow ~\~ • ■ • 



(4.10) 



One should note a striking similarity of the second term in Eq. ( |4.10i ) to the ^-term in Eq. (|4.9j). 



This gives some hope that the proposed approximate wave-functional (|4Jj) may in fact not only in- 
corporate Casimir scaling of string tensions at intermediate distances, but also Af-ality dependence 
in the asymptotic region. 

5. Conclusions 

• Casimir scaling is a natural outcome of a model of the QCD vacuum as a medium in which 
color magnetic fields fluctuate almost independently and are only weakly constrained by the 
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condition that the total flux through a cross-section of the domain conesponds to an element 
of the gauge group center. 

Lattice data show that string tensions of higher-representation potentials at intermediate dis- 
tances satisfy Casimir scaling to surprising accuracy for SU(2), SU(3), and even G2 gauge 
theory. 

Some elements of the model are (or may be) contained in the simple approximate form of the 
confining Yang-Mills vacuum wave-functional in (2+1) dimensions proposed in Ref. [21]. 
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